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Figure  1.  The  deployment  site  was  at  the  Field  Research  Facility  (FRF)  at  Duck,  NC. 


I.  INTRODUCTION 


An  array  of  five  acoustic  Doppler  velocimeters  (ADV),  which  produce  high  quality  measure¬ 
ments  of  the  three-dimensional  velocity  vector  in  a  sample  volume  with  a  scale  of  one  centime¬ 
ter,  was  deployed  from  late  August  through  late  November  of  1997  at  a  water  depth  of  ap¬ 
proximately  4.5  m  off  Duck,  North  Carolina  (Figure  1).  The  sensors  were  deployed  near  the  sea 
floor  but  above  the  centimeters-thick  wave  boundary  layer,  and  the  sampling  scheme  was  de¬ 
signed  to  resolve  turbulence  statistics  averaged  over  tens  of  minutes,  much  longer  than  typical 
wave  periods  but  shorter  than  time  scales  associated  with  variability  of  energetic  wind-driven 
and  wave-driven  alongshore  flows. 

A  part  of  the  SANDYDUCK  97  field  program,  the  experiment  addresses  turbulence  in  the  shal¬ 
low  nearshore  environment,  where  water  depths  are  on  the  order  of  meters,  energetic  currents  are 
forced  by  both  winds  and  breaking  waves,  and  the  motion  of  water  and  sand  is  of  primary  scien¬ 
tific  interest  and  practical  concern.  Existing  models  and  measurements  indicate  that  turbulence 
is  particularly  important  in  this  environment:  the  effect  of  bottom  drag,  which  plays  a  dominant 
role  in  controlling  the  magnitude  of  wind-driven  and  wave-driven  flows,  is  believed  to  be  trans¬ 
mitted  through  the  water  column  by  turbulent  Reynolds  stresses;  near-bottom  turbulent  stresses 
are  believed  to  control  the  entrainment  of  sediment  from  the  sea  floor;  and  turbulence  in  the  wa¬ 
ter  column  is  believed  to  maintain  sediment  in  suspension  against  the  action  of  gravity.  In  spite 
of  its  importance,  turbulence  in  the  nearshore  environment  has  rarely  been  measured.  Direct 
measurements  of  turbulent  Reynolds  shear  stress,  in  particular,  have  never  been  obtained  in  this 
environment,  because  of  fundamental  problems  produced  by  surface  waves. 

The  objectives  of  the  analysis  are  (1)  to  obtain  direct  estimates  of  turbulent  Reynolds  shear 
stress,  by  using  a  novel  method  involving  the  difference  between  velocity  measurements  ob¬ 
tained  by  pairs  of  spatially  separated  sensors;  (2)  to  obtain  indirect  inertial-range  estimates  of  the 
dissipation  rate  for  turbulent  kinetic  energy;  (3)  to  test,  in  collaboration  with  other  SANDY¬ 
DUCK  principal  investigators,  a  wave-averaged  alongshore  momentum  equation  in  which  wind 
stress  and  cross-shore  gradient  of  wave-induced  radiation  stress  balance  bottom  stress;  and  (4)  to 
test  a  simplified  turbulence  energy  balance  in  which  shear  production  balances  dissipation. 

Other  notable  components  of  the  SANDYDUCK  97  field  program,  for  the  present  purposes, 
were  an  array  of  pressure  and  velocity  sensors,  deployed  by  Steve  Elgar  (WHOI),  Britt  Rauben- 
heimer  (WHOI),  Tom  Herbers  (Naval  Postgraduate  School),  Bob  Guza  (Scripps  Institution  of 
Oceanography)  and  Bill  O’Reilly  (University  of  California  at  Berkeley);  an  Acoustic  Doppler 
Current  Profiler,  recording  flow  profiles  in  the  mid  and  upper  portions  of  the  water  column,  de¬ 
ployed  by  Peter  Howd  (University  of  South  Florida)  and  Kent  Hathaway  (FRF);  and  a  set  of 
wind  sensors  maintained  by  Chuck  Long  (FRF)  and  Jim  Edson  (WHOI)  on  a  mast  at  the  end  of 
the  FRF  pier  (Figure  2).  The  bottom-mounted  ADV  frame  was  deployed  between  the  two 
furthest-seaward  lines  of  Elgar,  Herbers,  Guza  and  O’Reilly,  so  that  an  estimate  of  the  cross- 
shore  gradient  of  wave-induced  radiation  stress  can  be  obtained  along  these  two  lines.  The  me¬ 
teorological  instrumentation  deployed  by  Long  and  Edson  provided  direct  covariance  estimates 
of  wind  stress,  approximately  20  m  above  the  water  surface,  throughout  the  measurement  period. 
A  NOAA  tide  gauge,  deployed  near  the  site,  provided  height  of  the  water  column  above  the 
frame. 
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The  purpose  of  this  report  is  to  document  the  instrumentation  and  deployment  of  the  ADV  array, 
to  present  an  overview  of  the  data,  to  summarize  the  data  processing  and  preliminary  analysis 
and  to  describe  the  formats  of  the  data  archives. 
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II.  INSTRUMENTATION 
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A  low-profile  frame  was  designed  and  built  to  provide  a  platform  on  which  to  mount  five 
SonTek1  acoustic  Doppler  velocity  (ADV)  meters.  The  instrument  array  (Figure  3)  included 
three  of  the  10  MHz  field  version  of  SonTek’s  acoustic  Doppler  velocimeter  (ADVF)  and  two  of 
the  more  rugged  5  MHz  ocean  version  (ADVO).  The  relative  positions  of  the  two  sensors  within 
each  sensor  pair  were  suitable  for  providing  estimates  of  Reynolds  stress  by  means  of  the  dif¬ 
ferencing  technique  (Trowbridge,  1998)  and  the  sensing  volumes  were  at  different  heights,  so 
that  the  measurements  can  provide  estimates  of  the  vertical  gradient  of  Reynolds-averaged  hori¬ 
zontal  velocity.  The  offshore  ADVF  (Figure  4)  was  likely  too  far  from  the  other  ADVFs  to  pro¬ 
vide  reliable  stress  estimates,  and  it  was  included  primarily  to  provide  information  about  vertical 
structure.  Measurements  from  either  the  pair  of  ADVOs  or  the  farthest-onshore  pair  of  ADVFs 
are  sufficient  to  achieve  the  objectives  described  in  the  previous  section. 


Figure  3. 
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The  upward-looking  ocean  probes  were  simultaneously  sampled  at  140  msec  intervals  and 
logged  on  a  TattleTale2  6-1M  for  approximately  26  minutes  at  the  beginning  of  each  hour.  The 
side-looking  field  probes,  ADVF_a,  ADVF_b  and  ADVF_c  were  logged  by  three  separate 
TattleTale  6F  loggers.  The  sampling  was  synchronized  by  using  a  master  logger  (ADVF_a), 
which  sent  a  sync  pulse  to  slave  #1  (ADVF_b)  and  slave  #2  (ADVF.c)  to  initiate  simultaneous 
sampling  and  logging.  The  field  probes  sampled  at  25Hz  for  9.6  minutes  at  the  beginning  of 
each  hour.  The  ADVs  were  placed  along  one  side  of  the  frame,  which  was  designed  to  face  the 
prevailing  along-shore  current  at  the  experiment  site.  The  battery  and  logger  pressure  cases  were 
placed  behind  ADVF_c,  where  the  sensing  volume  was  above  any  flow  obstruction.  The  heights 
of  each  sensing  volume  from  a  nominal  bottom  are  given  in  the  table  below. 


Sensor  ID 

Serial  Number 

Height  of  Sensing  Volume  (m) 

ADVF_b 

SN/1160 

0.36 

ADVF_a 

SN/1152 

0.46 

ADVO_B 

SN/5026 

0.76 

ADVCLA 

SN/0005 

0.86 

ADVF_c 

SN/1154 

0.98 

ADVO_A  was  equipped  with  a  strain-gauge  pressure  sensor,  compass,  tilt  meter  and  a  tempera¬ 
ture  sensor,  which  were  housed  in  the  probe  casing.  The  probes  were  coated  with  two  to  three 
coats  of  antifouling  paint,  except  the  transducer  faces,  which  were  coated  with  only  one  coat. 


A  static  salinity  and  temperature  were  defined  to  provide  identical  calibration  of  sound  speed 
(1529.7  m/s)  across  the  sensors.  Each  of  the  sensors  was  set  up  as  follows: 


Outformat: 

Temp: 

TempMode: 

VelocityRange: 

SampleRate: 

SampleMode: 

CoordiateSystem: 


BINARY 

20.31  Sal:  29.71 

USER 

ADVF  (250cm/s)  ADVO  (200  cm/s) 

ADVF  (25  Hz)  ADVO  (N/A) 

ADVF  (SYNC,  master  SAMPLE,  slaves)  ADVO  (SAMPLE) 
XYZ 


2 Onset  Computer  Corporation,  Pocasset,  MA  02559-3450 
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Figure  4.  Instrumentation  and  orientation 
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Figure  5.  Bathymetry  from  FRF  surveys  with  the  location  of  the  ADV  array  (X) 
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III.  DEPLOYMENT 


The  frame  was  deployed  on  8/25/97  off  the  coast  at  the  Field  Research  Facility  (FRF)  of  the  U.S. 
Army  Engineer  Waterways  Experiment  Station  (WES),  Coastal  Engineering  Research  Center 
(CERC),  using  the  Coastal  Research  Amphibious  Buggy  (cover  figure).  It  was  positioned  in  4.5 
meters  at  426.07  meters  offshore  and  830.38  meters  along-shore.  (See  Figure  2).  The  position 
was  noted  by  divers  and  is  relative  to  the  center  of  the  securing  post  which  was  nearest  ADVF_b. 
Three  separate  surveys  were  conducted  by  the  FRF  during  the  deployment  of  the  ADV  array, 
providing  the  bathymetry  near  the  deployment  site,  as  seen  in  Figure  5.  Five  posts  were  hydrau¬ 
lically  driven  approximately  15  feet  into  the  sandy  bottom  to  secure  the  frame,  both  above  and 
below  each  foot.  This  provided  a  stationary  height,  but,  as  noted  by  divers,  the  bottom  beneath 
the  frame  scoured  as  much  as  30  cm  (Figure  6a).  The  orientation  of  the  frame  was  observed  to 
be  such  that  the  instrumented  edge  of  the  tripod  was  80°  from  Magnetic  North.  The  magnetic 
deviation  at  the  site  is  approximately  10°  from  True  North,  which  confirmed  that  the  edge  was 
shore  normal,  as  planned.  Upon  recovery  of  the  instrument  (1 1/21/97),  it  was  found  that  the  base 
of  each  sensor,  the  channel  and  the  pressure  casings  were  encrusted  with  deposits  from  the  sand 
builder  worm,  sabellaria  vulgaris  (personal  communication,  A.  Frese  and  V.  Starczak  ,WHOI). 
(See  Figure  6b).  These  deposits  may  have  caused  significant  flow  disturbance  in  the  latter  part 
of  the  experiment. 

The  head  of  ADVO_A  developed  a  small  leak,  which  eventually  caused  intermittent  failure  of 
the  velocity  sensor  after  8/31  04:00  EST,  and  ultimate  failure  of  the  velocity  sensor  by  10/29 
16:00.  The  compass,  tilt  and  velocity  boards  were  affected,  but  the  pressure  appeared  to  be  pro¬ 
tected  and  continued  to  function  for  the  duration.  ADVF_b  was  bumped  on  8/31  04:00  which 
caused  contamination  of  the  vertical  axes  with  flow  from  the  x  direction  and  on  9/4  04:00  path  2 
dropped  out  (correlation  coefficient  went  from  about  96%  to  41%).  Upon  recovery,  the  angle  of 
tilt  seemed  to  be  on  the  order  of  20°  about  the  x  axis  with  the  shoreward  side  of  the  head  coming 
up  out  of  the  x-y  plane.  ADVF_a  was  completely  mangled  upon  recovery,  but  it  appeared  to 
work  well  through  9/3  19:00  EST.  The  field  probe  loggers  failed  after  10/10/97. 

Since  1978,  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)/National  Ocean 
Service  (NOS)  has  operated  a  primary  tide  station  (No.865-1370)  at  the  seaward  end  of  the  FRF 
pier,  tide  gauge  #11.  A  NOS  acoustic  tide  gauge  (Next  Generation  Water  Level  Measurement 
System,  NGWLMS)  is  used  to  collect  water  level  data  every  6  minutes  throughout  the  month. 
Data  are  referenced  from  meters  to  NGVD. 

Two  anemometers  were  mounted  on  a  mast  at  the  end  of  the  pier.  Jim  Edson  (WHOI)  deployed 
a  3-axis  Solent  ultrasonic  anemometer  model  1012R2A1  at  21.7  meters  above  NGVD,  providing 
wind  speed,  air  temperature,  tilt,  covariance  of  horizontal  wind  and  temperature  with  vertical 
wind  velocity.  The  data  were  observed  at  20.833Hz  for  twenty  minutes  continuously.  A  K-Gill 
Impellor  vane  anemometer  (Model  35301  )2  was  deployed  at  24.6  m  above  NGVD  by  Chuck 
Long  (FRF),  sampling  at  4  Hz  for  approximately  51  minutes  for  each  burst. 


-LSolent  mfgrd  by  Gill  Instruments  Limited,  Lymington,  Hampshire  (England) 
2K-Gill  mfgrd  by  R.M. Young  Co.,  Traverse  City,  MI 
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Figure  6b.  Fouling  of  frame  was  evident  upon  recovery  on  1 1/21/97. 


IV.  DATA  PROCESSING 


The  data  were  unpacked  and  each  burst  was  stored  in  a  Matlab®1  file,  using  the  month,  day,  hour  and 
minute  (EDT)  in  the  filename.  The  format  for  each  of  these  files  is  described  in  Section  VII.  The  data 
were  then  converted  to  real  world  coordinates  (Figure  4),  where  v  is  along-shore  southward  flow  (m/s), 
u  is  cross-shore  (shore-ward)  flow  (m/s)  and  w  is  vertical  flow  (upward)  (m/s). 

For  probes  ADVF_a,  ADVF_b  and  ADVO_B,  data  were  flagged  bad  when  the  correlation  coefficient 
was  less  than  70  and  were  replaced  with  NaNs.  Probe  c  had  signal  amplitudes  significantly  lower  than 
the  other  field  probes,  which  resulted  in  lower  correlation  coefficients.  Beam  3  of  field  probe  c  often 
had  correlation  coefficients  below  the  70%  threshold.  For  ADVF_c,  the  correlation  coefficients  were 
not  used  to  clean  up  the  burst  data,  but  they  were  used  in  determing  the  burst  averaged  statistics. 

Velocity  data  for  all  probes  were  wild  point  edited  by  replacing  all  points  greater  than  four  standard  de¬ 
viations  away  from  the  mean  of  each  burst  with  NaNs. 

For  ADVOJB,  about  91%  of  the  bursts  lost  no  more  than  5%  of  the  data  in  using  these  techniques. 
Most  were  flagged  bad  due  to  the  correlation  coefficient  test  and  these  were  primarily  between  year 
day  277  and  287,  when  the  signal  strength  was  low. 


The  correlation  coefficient  and  signal  strength  did  not  provide  a  reliable  means  to  determine  valid  data 
for  the  faulty  probe  ADVO_A,  so  a  technique  was  developed  which  compared  100  point  (14  second) 
blocks  of  ADVO_A  with  simultaneously  sampled  ADVO_B,  accepting  the  block  when  there  were  more 
than  10  valid  (non-zero)  samples  in  the  block  and  the  two  velocities  were  highly  correlated,  with  the 
squared  correlation  coefficient,  R~,  greater  than  0.9  and  0.8,  for  horizontal  and  vertical  velocities,  re¬ 
spectively.  Otherwise,  the  block  was  replaced  with  NaNs.  The  data  were  truncated  to  11000  points  per 
burst,  being  an  even  interval  of  100  blocks,  using  N  =  37:11037  and  38:11038  for  ADVO_A  and 
ADVO_B,  respectively,  to  provide  the  best  cross-correlation  between  the  probes. 

The  orientation  of  ADVF_c  was  resolved  by  determining  the  principal  axes  of  the  burst  averaged  ve¬ 
locities  during  the  first  six  days.  The  transformation  matrix  is  presented  below  and  was  applied  to  each 
burst  as  follows: 


[v  u  w]  =  (C'  x  [z  y  x]T 


c  = 


-0.5176 

-0.8437 

0.1423 


-0.8513 

0.4912 

-0.1845 


-0.0857 

0.2167 

0.9725 


1The  MathWorks,  Inc.,  Natick,  MA  01760 


9 


Two  sets  of  burst  statistics  were  created.  A  set  of  132  hours  of  data  between  8/25/97  10:00  and  8/30/97 
21:10  EST,  representing  the  mean  of  the  first  9.6  minutes  of  data  at  the  top  of  each  hour.  Times  when 
the  mean  correlation  coefficient  of  probe  ADVF_c,  path  3,  was  less  than  70  were  masked  with  NaNs. 
The  second  set  of  data  includes  25.7  minute  burst  statistics  from  the  ocean  probes  for  2083  hours.  Sta¬ 
tistics  using  data  from  ADVO_A  contains  NaNs  when  the  instrument  was  not  functioning  during  more 
than  90%  of  the  burst. 

The  standard  deviation  of  the  velocity  and  signal  strength  was  computed  from  detrended  data. 

The  mean  half-hourly  estimates  of  cross-shore  and  along-shore  bottom  stress  were  computed  using  a  dif¬ 
ferencing  technique  (Trowbridge,  1998).  Field  probes  were  paired  (ADVF_a  with  ADVF_b)  and  ocean 
probes  were  paired  (ADVO_A  with  ADVOJB).  Outliers  in  du,  dv  and  dw  (points  more  than  4  standard 
deviations  away  from  the  mean)  were  excluded  from  the  estimates  of  <u'w'>  and  cv'w'x 

As  documented  in  Section  VI,  our  pressure  sensor  did  not  provide  reliable  estimates  of  depth.  There¬ 
fore,  depth  was  derived  by  adding  the  tide  data  from  NOS  tide  gauge  #1 1  to  the  NGDV  elevation  of  the 
ADV  frame  (determined  by  the  location  and  the  FRF  surveys).  The  variance  from  pressure  was  not  af¬ 
fected  and  is  used  as  described  below.  Pressure  counts  were  converted  to  meters  using  Sontek  supplied 
calibration  coefficients:  meters  =  dbars  =  0.000287*counts-1.13. 

Dissipation  was  derived  from  spectra  computed  using  14  second  blocks  with  a  Hanning  window  and  no 
overlap.  Integration  was  performed  between  1  and  2  Hz,  in  an  attempt  to  stay  within  a  region  of  fre¬ 
quencies  above  the  wave  peak  and  below  the  noise  floor.  Spectra  were  corrected  for  advection  of  sur¬ 
face  waves  (Lumley  and  Terray,  1983;  Trowbridge  and  Elgar,  submitted). 

The  bottom  orbital  velocity,  U^,  was  computed  as  the  square  root  of  the  sum  of  the  variance  of  the 
horizontal  velocity  between  0.05  and  0.5  Hz.  Significant  wave  heights,  Hs,  were  estimated  by  convert¬ 
ing  the  pressure  spectra  to  surface  wave  variances,  St|T|,  using  linear  wave  theory,  and  integrating  in  the 
wave  peak  (0.05  -  0.21  Hz).  Significant  wave  heights  were  also  calculated  from  horizontal  velocity  vari¬ 
ance,  using  linear  wave  theory,  and  are  compared  with  the  estimates  derived  from  the  pressure  data  in 
Section  VI.  For  these  parameters,  spectra  were  computed  using  1028  point  (144  second)  blocks  with  a 
Hanning  window  and  no  overlap.  The  mean  frequency  in  the  wave  band  was  computed  as  a  weighted 
mean: 

fa5Hz  /  \ 

Jo.o5Hz(Suu+Svv)>f<df 

^bar  /.0.5Hz 

Jo.05Hz(SUU  +  SVV)  *  df 

The  sonic  anemometer  data  were  converted  to  the  along-shore,  oceanographic  coordinate  system  (Fig¬ 
ure  4.).  The  twenty  minute  burst  statistics  were  computed  without  detrending  or  windowing.  Variance 
greater  than  0.6  m2/s2  was  replaced  by  NaN.  The  time  (year  day)  was  converted  to  EST,  with  0.5  desig¬ 
nating  noon  on  1/1/1997  and  linearly  interpolated  to  correspond  with  the  9.6  and  25.7  minute  burst  data 

for  each  set.  Density  (p)  was  computed  using  air  temperature  (T),  p  =  101300/287/(T+273.16)  and  heat 
flux  was  derived  as  1004*p*<wT7>. 
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V.  DATA  SUMMARIES 


Figures  7-19  present  data  from  all  five  ADVs  along  with  wind  statistics  from  the  sonic  anemometer, 
for  the  first  five  and  a  half  days,  when  all  the  sensors  were  functioning  well.  Figures  20  -  27  present 
data,  as  it  was  available,  through  the  first  60  days.  The  data  were  truncated  for  the  report,  since  there 
were  no  turbulence  statistics  available  after  ADVO_A  failed. 

The  significant  wave  heights  presented  are  those  computed  from  the  pressure  spectra  and  the  bottom  or¬ 
bital  velocity  are  those  computed  from  horizontal  velocities.  For  the  long  time  series,  the  bottom  orbital 
velocities  presented  are  those  from  the  horizontal  velocities  of  ADVO_B. 


11 


Figure  7.  Unfiltered  flow  (up  is  along-shore,  southward  &  right  is  onshore) 
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Figure  9.  Alongshore  velocity 


Figure  12.  Bottom  orbital  velocity 
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Figure  20.  ADVO_B  velocity 
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Figure  23.  Estimates  of  dissipation  (1 0  W/kg) 
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Figure  26.  Wind  speed  from  sonic  anemometer 
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VI.  DATA  ANALYSES 


Comparison  of  wave  characteristics  from  ocean  probes  and  pressure  sensor 

Figure  28a  shows  that  the  significant  wave  heights  computed  from  velocity  and  pressure  are 
highly  correlated.  Figure  28b  documents  the  strong  correlation  in  the  wave  peak  between  the 
two  ocean  probes.  The  data  are  only  compared  for  the  first  5.5  days  when  both  probes  were  fully 
functional  and  the  full  wave  peak  could  be  resolved. 


Figure  28a.  8/25/97  -  1 0/23/97  Figure  28b.  8/25/97  -  8/31/97 


H  from  pressure  (m)  H  from  ADVO_A  (m) 


Comparison  of  depth  measurements  from  FRF  gauge  and  ADVO_A  pressure  sensor 


Comparison  of  pressure  from  ADVO_A  with  the  FRF  tide  gauge  No.  11 


Tide  gauge  data  from  NOS  Tide  gauge  11,  averaged  to  represent  the  first  half-hour  of  each  hour, 
are  compared  here  with  pressure  data  from  the  strain  gauge  sensor  mounted  in  the  end  cap  of 
ADVO_A.  Figure  29  documents  the  offset  presumably  from  the  temperature  sensitivity  of  the 
strain  gauge  sensor.  Figure  30  compares  the  tidal  fluctuations  and  indicates  that  the  temperature 
sensitivity  of  the  strain-gauge  pressure  sensor  only  affects  an  offset,  not  a  change  in  gain.  Low- 
pass  filtering  was  applied  using  p!64t.m  (Rosenfield,  1983)  to  derive  the  tidally  averaged  depth, 
so  it  could  be  removed. 


Figure  30.  Tidal  variance  in  depth  from  FRF  gauge  vs.  ADVC^  pressure  sensor 


depthtgii  “  Pl64t(depthtg)  (m) 
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Comparison  of  wind  data  from  FRF  K-vane  anemometer  with  Edson  sonic  anemometer 

The  data  from  the  sonic  anemometer  are  compared  here  with  the  data  from  the  vane  anemometer. 
Much  of  the  scatter  in  the  comparisons  of  the  turbulence  statistics  can  be  attributed  to  different 
processing  techniques.  The  vane  data  were  processed  by  detrending  the  data,  using  smaller  win¬ 
dows  (512  second)  and  cleaning  up  times  when  there  was  slow  wind  speed  or  the  wind,  at  any 
time,  passed  through  the  5-degree  ’dead  zone’.  As  seen  in  Figure  32,  both  sensors  observe  en¬ 
hanced  turbulence  when  any  component  of  the  wind  blows  off-shore. 


Figure  32.  Comparison  of  observed  stress  with  estimated  stress  (Large  &  Pond,  1981) 

Cross-shore  wind  stress  (Pa)  Along-shore  wind  stress  (Pa) 
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VII.  FILE  DESCRIPTIONS 


ADV  FIELD  PROBE 
Binary  data  flies 

All  three  loggers  placed  a  time  stamp  at  the  beginning  of  each  burst.  Each  file  block  is  529288 
bytes  long  and  contains  two  bursts.  Each  burst  contains  data  which  were  logged  continously  for 
9.6  minutes  at  25  Hz  and  formatted  as  described  below. 


ADV  field  probe  format 

Each  block  began  with  a  date  stamp  (EDT):  MMDDHHMNSE  (month, day,hour, minute, second) 
and  was  followed  by  two  sets  of  about  14400  records  as  follows: 

Variable 

#  variables 

bytes/variable 

bytes 

units 

Keyword(81 12) 

1 

2 

2 

Sample  Id 

1 

2 

2 

1  -  number  of  samples  in  burst 

Velocity 

3 

2 

6 

(0.1  mm/sec) 

Signal  Strength 

3 

1 

3 

(counts,  where  0.43  dB/count) 

Correlation  Coefficient 

3 

1 

3 

(0  to  100,  where  >  70  is  considered  ok) 

Checksum 

1 

2 

2 

sum  of  bytes  plus  base  (0xa596) 

Total  bytes/record: 

18 

The  data  were  unpacked,  and  truncated  to  the  first  14400  samples,  using  binadvfps.m  and 
docleans.m,  where  s  is  the  sensor  id  (a,  b  or  c). 

Processed  burst  data  files 

Field  probe  data  were  stored  on  one  CD  for  8/25  11:00  -  9/4  23:00  (EDT).  Each  sensor’s  data 
was  stored  in  its  respective  file  which  was  named  sMODAHRx.mat,  where  s  is  a,  b  or  c  for 
ADVF_a,  _b  &  _c,  respectively,  and  x  is  either  V  for  velocity  data,  A  for  signal  strength  or  C  for 
correlation  coefficients.  These  files  include  the  following  data: 

vs,  us,  ws  -  along-shore,  cross-shore  and  vertical  velocity  (m/s), 

where  s  is  the  sensor.  Data  includes  NaNs. 
samp  -  signal  strength  (dB)  of  sensor  s  (a,  b  or  c) 
see  -  correlation  coefficient  for  each  of  the  sensors  (s)  (%) 
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ADV  OCEAN  PROBE 


Raw  data  files 

The  ocean  probes  were  logged  simultaneously.  Each  140  milliseconds,  a  burst  was  recorded 
from  each  sensor,  beginning  on  the  hour  and  recording  for  25.6  minutes.  Two  bursts  were  re¬ 
corded  per  file  block  (1015808  bytes). 


r  ADV  ocean  probes  format 

Each  block  began  with  a  date  stamp  (EDT):  MMDDHHMNSE  (month,day, hour, minute, second) 
and  was  followed  by  1 1038  records  as  follows: 

Variable 

#  variables 

bytes/variable 

bytes 

Probe 

units 

Key  word(  871c) 

Sample  Id 

Velocity 

Signal  Strength 

Correlation  Coefficient 

heading 

pitch 

roll 

temperature 

pressure 

Checksum 

Keyword(8 1 12) 

Sample  Id 

Velocity 

Signal  Strength 
Correlation  Coefficient 
Checksum 

1 

1 

3 

3 

3 

1 

1 

1 

1 

1 

1 

1 

1 

3 

3 

3 

1 

2 

2 

2 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

1 

1 

2 

2 

2 

6 

3 

3 

2 

2 

2 

2 

2 

2 

2  : 
2 

6 

3 

3 

2 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

B 

B 

B 

B 

B 

B 

1  -  number  of  samples  in  burst 
(0.1  mm/sec) 
(counts,  where  0.43  dB/count) 
(0  to  100,  where  >  70  is  considered  ok) 
(0.1  degrees) 
(0.1  degrees) 
(0.1  degrees) 
(0.01  degrees  C) 
(counts) 

sum  of  bytes  plus  base  (0xa596) 

1  -  number  of  samples  in  burst 
(0.1  mm/sec) 
(counts,  where  0.43  dB/count) 
(0  to  100,  where  >  70  is  considered  ok) 
sum  of  bytes  plus  base  (0xa596) 

Total  bytes/record: 

46 

Processed  burst  data  files 

These  unpacked  raw  files  were  cleaned  up  (docleanA.m  &  docleanB.m)  and  shifted  so  that 
Ocean  Probe  A  &  Ocean  Probe  B  are  in  phase.  (See  Section  III.)  The  cleaned  up  data  files  con¬ 
tain  vP,  uP,  and  wP,  where  P  is  either  A  or  B  in  files  named  sdP_MMDDHHV .mat,  and  the  cor¬ 
relation  coefficient  and  signal  strength  (in  integer  dB)  were  stored  in  sdP_MMDDHHC  or 
sdP_MMDDHHA.mat,  respectively.  The  data  for  Ocean  Probe  A  should  be  dealt  with  in  blocks 
of  100,  since  the  processing  was  completed  on  those  blocks.  (See  dissipation.m.)  The  pressure 
(m)  and  temperature  (degrees  C)  are  stored  in  sdP_MMDDHHPT.mat.  Time  (month,  day  and 
hour)  for  these  files  are  taken  from  the  filename  (EDT)  and  the  minute  can  be  reconstructed  by 
adding  (37*0. 14:0. 14: 11000)760/24  to  the  day  of  year. 
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BURST  AVERAGED  PROCESSED  DATA 


By  computing  the  burst  averages  of  each  parameter,  hourly  estimates  centered  at  12.8  minutes 
after  the  top  of  the  hour  are  stored  in  sandyduckL.mat.  Time  (tdoy)  is  year  day,  with  0.5  repre¬ 
senting  noon  on  January  1,  1997  (EST).  windL.mat  contains  the  sonic  anemometer  hourly  aver¬ 
aged  wind  statistics  for  the  long  record. 

sandyduckS.mat  contains  the  first  six  days  only,  containing  the  field  probe  data  and  the  burst 
averaged  data  from  the  ocean  probes,  during  the  times  corresponding  to  when  the  field  probes 
were  sampling.  Time  (tdoy)  is  year  day  (EST),  with  0.5  representing  noon  on  January  1,  1997. 
windS.mat  contains  the  wind  statistics  from  the  sonic  anemometer,  interpolated  to  correspond  to 
the  burst  averaged  time  of  the  field  probes. 
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